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Abstract 

The relation between two versions of so called non-equilibrium statistical operator 
method (NESOM), NESOM-1 due to Zubarev (1961) and NESOM-2 due to Zubarev 
and Kalashnikov (1970), is considered. It is proved that, once the balance equations 
of NESOM-2 are satisfied, those of NESOM-1 will be satisfied with the same set 
of the macro-parameters. The proof uses the convexity-type inequalities, and does 
not involve any assumptions additional to the rationales behind NESOM. However, 
converse statement cannot be proved within this technique. An extension of the 
proof to overall equivalence of the two nonequilibrium ensembles is discussed. 
Dedicated to the memory of Dmitrii Nikolaevich Zubarev 



Preface 



The methods of non-equihbrium statistical ensembles based on maximum entropy (Max- 
Ent for short) principle are nowadays well documented. One of them, developed inde- 
pendently by Zubarev [1] and McLennan [2] extended the Gibbs ensembles method to 
non-equilibrium. Later method, due to Zubarev and Kalashnikov [3, main text], [4, main 
text] made also use of the motion quasi-invariant construction [1] but applied it to the 
MaxEnt distribution rather than to the entropy itself. Both developments are referred to 
as Zubarev's Nonequihbrium Statistical Operator Method (NESOM), while the notations 
NESOM-1 for the earher and NESOM-2 for the later version of NESOM, respectively, 
were introduced to distinct them [6]. Other methods, which are close to NESOM, were 
pubhshed in the interim [3] , [4] , [5] . 

When NESOM-2 became a method of choice, the expected question of its equivalence 
to NESOM-1 was promptly addressed that time [6]. Other proofs of the equivalence, 
restricted to classical systems, were published later [7], [8]. More later, the present author 
and Kalashnikov showed that the above mentioned proofs used assumptions, not found 
among the NESOM rationales, and proposed another proof [9] , which seems satisfactory 
from the viewpoint of rigor adopted in physics. Though available in Enghsh, Ref.[9] 
remained un-noted - the recent review [10] refers the proofs of Refs.[6], [8] to as the 
ultimate results on the problem. For that, and not only that, reason I decided to republish 
Ref.[9] on Internet. My motivation for this is better explained by the passage below. 

Albeit much more applications were treated using NESOM-1 than NESOM-2, no work- 
ing perturbation technique beyond lowest-order approximation (such as high-temperature 
or weak-interaction approximations) were developed in the former method. Moreover, be- 
cause of rather involved statistical distribution structure in NESOM-1, such a development 
seems desperate. On the other hand, the statistical distribution construction in NESOM-2 
is certainly analogous to the scattered wave function construction in the Gell-Mann and 
Goldberger form of scattering theory, the MaxEnt distribution being an analog of the in- 
coming wave function. Due that analogy, the Gell-Mann-Goldberger integral equation based 
perturbation technique applies to NESOM-2 [3, main text]. So, once the equivalence of 
NESOM-1 and NESOM-2 is proved, there is no more need to deal with cumbersome and 
intractable perturbation series of NESOM-1. 

Just within the NESOM-2 framework, on such important examples as plasma screen- 
ing and Kondo effect at non- equilibrium, partial summation of the perturbation series was 
shown to be performable [11]. However, this trend was not focused on since then. As 
a result, in spite of its conceptual advantages {see review by Luzzi et al. [10]), NESOM 
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legs behind of regular methods {such as the Keldysh and Kadanoff-Baym ones) regarding 
the state-of-art respect. Nevertheless I believe that the potential of NESOM-2 {contrary to 
NESOM-l) has not yet been exhausted and that a proper diagrammatic technique devel- 
opment would benefit NESOM. 

The present publication is rather a 'remake' than reproduction of Ref.[9]. The text of 
the paper is essentially revised but the references remained in the format of the original. 
(In Preface the citation of Russian papers and books, even cited in the paper text, done 
on their English translations, when available). Upon revision I did not aimed to make 
the proof as rigorous as the modern Mathematical Physics standard requires [12], since 
nowadays the level of rigor customary for Equilibrium Statistical Mechanics [13] still 
remains unreachable for most problems of Nonequilibrium Statistical Mechanics. 
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1 Introduction 



NESOM is formally based on special construction, which involves the so called quasi- 
equilibrium statistical operator (QESO) 

p, {t, 0) = e-^(*'°) , 5 (i, 0) = $ (t) + J2Fn (t) Pn . (1) 

n 

Here t is the time variable (which is dummy one unless dynamics enters the play), {Pn} 
is a set of gross variables, that is observables expectation values of which describe the 
non-equilibrium state of interest at instant t. $ (t) is the Massieur-Planck function 



^{t) =lnTr 



-EF„(t)Pn 

e " 



(2) 



and Fn {t) are the macro-parameters conjugated, in thermodynamical sense, to the gross 
variables averages {Pn)l ■ The first set of the equations, connecting F„ (t) and {PnYg, 
follows from Eq.(2) 

(Pn)l^Tr[p,{t,0)Pn]^-^^y (3) 

To arrive at the second set, define over the whole set of density matrices p the information 
entropy functional 

e[p]^-Tr{p\np) , (4) 

and the quasi-equilibrium entropy 

E (i) = € [p,] = {S {t, 0))l = $ (t) + 5^ F„ (t) {Pn)l ; (5) 

n 

thus S {t, 0) may be called the entropy operator. Then it is proven that 

Fn{t) = ^. (6) 

Consider now the dependence of the operators on time. Any operator Q may explicitly 
depend on t that is notated by Q {t, 0) (this notation has already been used above for 
QESO and the entropy operator). Operator Q (t, 0) is said to be dynamics invariant (or 
integral) if it satisfies the Liouville equation 

Q it, 0) = + i£ it) Qit, 0) = , (7) 
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where £ (t) is the Liouville super-operator, which action is given by the commutator (in 
classical case Poisson-bracket) with the Hamiltonian. In particular, any non-equilibrium 
statistical operator (NESO) satisfies Eq.(7) and QESO does not. One of the crucial 
concepts in NESOM is the dynamics quasi-invariant (or quasi-integral) . Given an operator 
Q {t, 0) and a positive number e, define the dynamics quasi-integral associated with Q {t, 0) 

by 



t t 



Q{t^) = e J e^(*°-*)Q {to, to - t) dto ^Q{t,0)- J e^^*""*) Q {to, to - t) dto , (8) 

— oo — oo 

where Q{ti,t\ — t-i) — il(ii, ^2) Q (^i, 0) and 11(^1,^2) is the dynamic evolution super- 
operator. It satisfies the couple of dual to each other equations 

^il(ti,t2) =11(^1,^2) iii(ii) (9) 



and 



J-il(^l,^2) = -^ii(^2)ii(^l,^2) (lo) 
OT2 



with the initial conditions il(t2,^2) = il(^ij^i) = ^, where and is the identity super- 
operator. For the time-independent Hamiltonians, Eqs.(9), (10) give il(ii, ^2) = 6**^*^"*^^''^. 

£ 

As seen from Eqs.(7) and (8), Q{t,Qi) — Q{t,Qi) if Q{t,Qi) is the dynamics invariant. 
Moreover, Q {t, 0) satisfies the equation 

d 



^ Q (t, 0) + z£ {t) Q (t, 0) = s 



g(t,o)-g(t,o) 



11) 



approaching to Eq.(7) at £ J, that explicates using the term 'quasi-invariant'. 

For a system occupying a finite domain Q of volume \ NESOM builds from QESO so 
called quasi-non-equihbrium statistical operator (QNESO) using the procedure of Eq.(8). 
QNESO is denoted here by the notation ("^p^ {t, 0), where the superscript a = 1,2 points 
to the NESOM version. More historically earlier one (NESOM-1) was formulated in the 
papers of D.N. Zubarev [1] and McLennan [2] as the generahzation of canonical ensemble 
method to non-equilibrium. In terms of Eq.(8) QNESO-1 is quasi-canonical distribution 



of the form 



(1), 



)p, (t, 0) = e*^W-^(*'°) , (t) = - InTr 



-5(t,0) 
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The later version (NESOM-2) was developed by D.N. Zubarev and V.P. Kalashnikov [3]. 
QNESO-2 is the dynamics quasi-integral built from pq(t, 0), that is 



(%(t,0)=p,(t,0) = e-^(*.o). 



(13) 



Eq.(12) includes additional normalization to the unity trace. The logarithm of the 
QNESO-1 trace normalization factor ^£(t) plays a crucial role in the present paper. In 
Eq.(13) the normalization holds automatically. 

The equation for ^'^^p^ {t, 0) results directly from Eq.(ll): 

d 



dt 



+ i£ (t) 



(2), 



W{t,0)=e[pg{t,0)- %e{t,0)] 



(14) 



More involved equation for ^-^^p^ (t, 0) is obtained by applying to the operator exponent in 
Eq.(12) the rules of time differentiation and of £ {t) [4], and then Eq.(ll). The result is 



^% {t, 0) 
= ej Wp,(i,0)e^^'(*^ <! Wa! 



S{t,0) -S{t,0) 



(15) 



where 



by the definition and 



(°)(^)* ^Tr [(")p,(t,0) A] 



(16) 



(17) 



The source-like terms emerged in the rhs of Eqs.(15), (14) break time-reversal invariance 
of the Liouville equation. Non-equilibrium may be described using expectation values 
of the type given by Eq.(17) while e tends to zero. But, due the Bogoljubov's idea of 
quasi- averages, any invariance-breaking perturbation is switched off only after performing 
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the thermodynamic limit (TL). So the quasi- averages, which would sustain the broken 
symmetry of time arrow towards future, may be defined by 



(a) 



{Ay = Urn Urn ^'^\A)1 . (18) 



Of course, Eq.(18) makes any sense if A is an intensive observable (hke density of an 
extensive observable) that may be assumed without any loss of generality. As was argued 
by D.N. Zubarev [4], such a subordination of the e limit (eL) and TL is crucial for the 
correct simulation of irreversibility.^ 

Eqs.(3) and (6) by no means define the macro-parameters. The ultimate rationale of 
NESOM, which facilitates a closed description of non-equilibrium, is the balance equations 
(BE). BE may be presented in different equivalent forms. The basic one is the self- 
consistency form that reads 

^^^HPn)' = (Pn)' ^ lim (P^, (19) 
|S2|— »oo 

but actually the differential form of BE 

|w*=(")^A(t,o)y (20) 

is exploited in the practice. Shghtly different route consists of using the 'pre-hmit' BE 

^''HPn)l-{Pn)l (21) 

and 

|w: = ^"^(A(t,0)y , (22) 

at yet finite but large |Q| and small e, and then performing TL and eL in due order (see 

D.N. Zubarev's review in Ref.[4]). For a = 2, Eq.(20) and Eq.(22) exactly follow from 

Eq.(19) and Eq.(21), respectively, (see proof in e.g. [4]) while for a = 1 the question 

of relation between the two forms of BE is not so simple. A discussion on this issue is 

presented in Appendix 1. 

^TL and eL may be performed simultaneously, but anyhow £ |r2| ^ oo [1], [1, Preface] 



As seen, BE are complicated (non-local) functional equations for defining F„ (t). NE- 
SOM asserts that solutions to BE, *^"'F„ (t), do exist and finds them in physical approxima- 
tions. In the first order with respect to the entropy production, QNESO-1 and QNESO-2 
are the same [4]. Thus, NESOM-1 and NESOM-2 are equivalent in that approximation. 
This concerns both the coincidence of the BE solutions 

(^)F, it) = it) (23) 

and of the averages over the two NESOM ensembles 

(24) 

at least for some class of intensive observables. While Eq.(24) holds, Eq.(23) holds au- 
tomatically, but not vice versa. Quite naturally, the general equivalence of NESOM 
ensembles is in question. Authors of Ref.[5] discussed a condition for Eq.(23) to hold. 
In Refs. [6] and [7, Preface] proofs of Eq. (24) were suggested. As shown in Appendix II, 
these proofs essentially used assumptions which, in fact, are not found among the ba- 
sics behind NESOM. Moreover, it is shown that the condition of Ref.[6] and 'asymptotic 
normalization' condition of Ref . [7, Preface] can not be satisfied. 

In this paper a proof of Eq.(23) is presented. This proof is based only on the NESOM 
rationales and on a number of assumptions, which seem to be tacitly adopted in NESOM 
and in other methods of reduced description of non-equihbrium (see e.g. [3, Preface] - [5, 
Preface] ) . 

2 Pierls-Bogoljubov Inequality 

The main tool in this paper is the Pierls-Bogoljubov inequality [8] (also [13, Preface]). 
This inequality is used in the two-sided form: 

Tr\e-^(A-B)] , , Tr [e"^ (A - 5)1 , , 

; — ^ < In Tr (e-^) - In Tr (e"^) < — , , 25 

Tr e-^ - V ; \ J- Tr e-^ ' ^ ^ 
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the equality being possible only a,t A — B. To begin with, put A — — In [("^p^ (t, 0)] and 
B = S {t,0). On account of the trace normahsation, the left side of Eq.(25) leads to the 
inequality 

(")(5(t,0))* >e[('^k] ^ ('^)e,(t) (26) 
that, upon satisfying BE in the form of Eq.(21), transforms to 

S(t)> (")6,(t), (27) 

where the Ihs also becomes dependent on e (and generally on a). This inequality expresses 
the MaxEnt principle at finite \Vt\ and e. The function ^^^teg (t) may be called quasi-Gibbs 
entropy. Post TL and eL form of the MaxEnt principle, the only feasible for the NESOM's 
route with Eq.(19), holds for the specific entropies 

a (t) = (t) = hm ("Ve (t) > hm ("^5, (t) = ("^5 (t) , (28) 



where 



are the specific quasi-equilibrium and pre-£L non-equilibrium entropies, and 

("^6 (t) 

(").,(t)= lim (30) 

is the specific pre-sL quasi-Gibbs entropy. It is assumed that all limits in Eqs.(28) - (30) 
exist and BE are satisfied (see next section). 

3 Handling TL and eh 

TL and eh are inevitable ingredients of NESOM. In order to accurately resolve TL issue 
(but, at the same time, do not enter its subtleties [8], [13, Preface]) it is expedient to 
introduce several relevant axioms, which are necessary for the NESOM formulation to 
make any sense and seem quite satisfactory from the viewpoint of rigor adopted by physi- 
cists. As to eL, it is treated here in a simplified manner (for a mathematically correct 
treatment, see [7, soon on Web]). 
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Axiom 1 The index n is discrete at finite \Q\ and continuous at \Q\ oo, while there 
exists a measure such that 

|ii^^pE^'^ = / Snd^r^ (31) 

for every hounded function Qn- The Fn{t) by {PnY and by lim ^"\Pn)\ products are 
^-integrable? 

Axiom 2 $ {t) is an extensive quasi-thermodynamic potential, in the sense that there 
exists the limit 

|f2|— »oo 

which is an analog of the equihbrium specific grand-canonical potential [9]. The above 
axioms assure existence of the functions defined in Eqs.(28), (29). Firstly ^{t), being an 
analog of the equilibrium specific entropy [9] , is given by 

a{t)^(P (t) + jFn (t) {PnY d^in . (33) 

Post TL form of Eqs.(3) and (6) express (i-*n)* and F„ (t) by variational derivation of (t) 
and a {t) in F„ (t) and (-Pn)*, respectively (keeping t frozen), so that (j){t) and a (t) are 
connected to each other via Legendre transformation, e.g. 



<l^{t)^<t)- J (P„)*^^(i/.„. (34) 
Lastly ("Ve (t) and ^"V (t) are given by 

(")(7, (t) ^cl>{t)+ [f^ (t) hm (")(P„)* dpin (35) 



and 

("V (t) ^Ct>{t)+ j Fr, (t) (P„)* dfXn , (36) 

respectively. 

^Although in n there may be a continous species before and remain a discrete species after TL, this 
axiom may be assumed without any loss of generality. 
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Axiom 3 lim (P„ (to, to — ^i))*^ exists and (by physical reasons) is continuous bounded 

|r2|— >oo 

function of time arguments.^ The product Fn (to) by lim {Pn(to,to — ti))*^ is fi-integrable 

|n|— >oo ^ 

function of n. 

The following corollary ensues from these axioms and the rationales of NESOM. 



Corollary 1 lim (^^(P„ (to, to — t))l exists and the product Fn (to) by '''^\Pn (to, to — t))l 
>oo 

is /i-integrable function. In addition 

^'\P„(to,to-t))'^ (^)(P„)*° . (37) 

Proof. By the definition of Eq.(17) and the evolution super-operator property 

^'\Pn(to,to-t))l^Tr{Pn[iX(t,to) ^(t,0)]}. (38) 

The action of the super-operator e^^*~*°^iL(t,to) on both sides of Eq.(14) and the integra- 
tion of the resulting equation from to to t leads to 

t 

e'^'-'°^iX(t,to) %e(t,0)= ^(to,0) + e j e'^''-'°^ii(t,,to) p,(ti,0) dti . 

to 

The use of this operator identity and Eq.(38) gives the following functional identity 

t 

(2)(P„ (to, to - t))l = (2)(P„)? e^(*°-*) +eJ{Pn (to, to - ^i))*^ e^^'^'^Ut^, (39) 

to 

from which the lemma statement follows, including Eq.(37). Note that performing eL in 
Eqs.(35) and (39) is quite rigorous within the NESOM rationales. ■ 

Though assuming existence of ^^^^Se (t) is physically reasonable, nothing assures this 
formally. The statements below connect that question for a = 1 to the behavior of ^'^ (t) 
in TL. 



^E.g. for bounded (in the Hilbert space of the system) operators P„, which norms \\Pn\\ are independent 
of Vq, the absolute value of this function is bounded by ||-Pn||. 
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Lemma 2 (a) Existence of^^^5e{t) is equivalent to extensivity of'^^it), that is to exis- 
tence of 

(b) For both '4>e (t) and ^^''Se (t) to exist it is sufficient that lim ^^^(Pn {to, to — t))l exist 

\Q\^OD 

and the F„ (to) by lim ''^\Pn (to, to — t))* products are fi-integrable. 

|Q|-+oo 

Proof. Differentiation of the definition of '^^ (t) in Eq.(12) gives the relation 

W{S(i,0))' - W/SM) \ . (41) 



By the definition of Eq.(26) and Eq.(41) 

(t) = (t) + (^S^) ^ = S {t, 0))l - (^1 + e-'^^ (t) . (42) 

This equation may be integrated to give 

t 

^e{t)^e J e^(*°-*) p\S{toM° - ^'^&eito)] dto. (43) 

— oo 

Due to Eqs.(26), (43) (t) > 0, the fact that will further be proved in other way. Eqs.(42) 

and (43), due to extensivity of (5* (t, 0))* expressed by Eq.(35), show that existence of 
-06 (t) does provide that of ^^^s^ (t) and vice versa. Eq.(8) and Axiom 1 give 



1(1) ' ^ ^ * 



J ^e(t,-t) j (1) (^^^ _ ^)^t ^^^^^^^ (44) 



— oo 



provided that the condition of (b) holds. As seen from the above consideration, Eq.(44) 
assures existence of both ip^ (t) and ^-^^s^ {t). On physical level of rigor, the condition of (b) 
may be adopted as being also necessary for (a). It is worth emphasizing that an attempt 
to prove the analog of Corollary 1 for ^^\Pn (to, to — t))l, making use of Axioms 1-3 alone, 
fails. ■ 

In TL Eq.(43) and the first equality in Eq.(42) transform to 

t 

i^e (^) = ^ / [^'Vs {to) - ^'he {to)] dto (45) 
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and 

t t 
^%{t) = s J e^(*°-*)a(io)c^to-^.(t) + J J F^{to)^p^dii^dto, (46) 

—oo — oo 

respectively. Unfortunately nothing may be inferred from Eqs.(45), (46) on existence of 
(^)s (t) and 

V'(i)=limV'.(^), (47) 

£j,0 

except the fact that these quantities, if exist, should be independent of time: ip {t) = ip 
and (1)5 (t) = (1)5. 

4 Equivalence of Two NESOM Ensembles 
4.1 Generating Functional in NESOM-1 

Let f = {/„} is a set of functions independent of time, satisfying the same requirements 
as Fn {t) . Define upon f the functional 

*e {t; f) = - In Tr [e-^^(*'f)] , (48) 

where 

S, {t; f) = S^L^f) + fnPn . (49) 

n 

Formally, (t; f) is the logarithm of trace normahsation factor for the statistical distri- 
bution with an 'entropy' defined by Eq.(48). While the use made of the last equality in 
Eq.(8), this 'entropy' is presented as 



t 

{t; f) {t, 0) + ^ fnPn - J e^^*°-*) S {to, to - t) 



dtn 



that may be thought as an artificial perturbation, by the shifts F„ (t) — > F„ (t) + fn, 
concerning the QESO entropy part with no effect on the entropy production. The distri- 
bution under consideration is auxiliary one, which tends to 'shifted' QESO, while zeroing 

the entropy production. The Taylor series for [t; f) has the form 

00 ^ I 

{t- f) = w + E /! E ^i'^ -' n (50) 

;=1 n\...ni i=l 
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where (t) was defined earlier and 



d'^e {t;f) 



dfn^-dfm 



I > 1. 



f=0 



For particular I — 1, 



(51) 



(52) 



Higher derivatives of (i; f) at f = are connected with the cumulant correlators of 
gross variables over QNESO-1 



(53) 



where 
« {A, . .. 
with 



1 1 i-i 



(54) 







(55) 



and the 'irreducible' ordered products calculated via the following recursive relations 



[A, (AO: ...Am{Xm)l]l=[A, (AO: ...Am{Xm)T-^'H[A, iX^)l ...A^iX^)l] 



m— 1 



-EE (^n): - A. (A,j:],^^ ( {K,^X ...A. (A, 



fc=l Trr 



(56) 



In these relation tt^ are all permutations of the type 



ll...lk lk+l---lm 

l....k k + l...m 



; ii < ... < ik, ik+i < ■■■ < ir. 



(57) 



and [B^ (Ai) ...Pp(Ap)]^ is the Dyson-ordered product, in which the factors with larger 
A's are stood to the left. Note that (^4 • B)\ is a non-equihbrium analog of the Kubo- 
Duhamel correlator [8] (also [13, Preface]). In spite of apparent assymetry relative to Ai 
of {A^ ■ ... ■ Ai^i ■ Ai)\, as defined by Eq.(54), this correlator appears to be symmetric 
with respect to all permutations of its operator constituents. 
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Consider the limit 



(t; f) = lim (5^ 

|S2|— >oo Si 



Term-by-term analysis of the series in Eq.(50) based upon Axiom 1, Eqs.(52) and (53) 
leads to the conclusion that (t; f) would be represented by functional Taylor series 



{t; f)=^e{t)+f2^ f ... f (t; ni, m) fl fn,d^n, , 
1=1 -I i=l 



(59) 



where existence of the first and second terms is provided by adopting Eq.(40) and by the 
rationale of NESOM-1, respectively. The higher-order terms would be made meaningful 
by assuming existence of limits 

« [P„, • ... • Pj^ = lim (P„, • ... • P„J* , (60) 



SO that 



and 



xW(t;n)= hm ^ (P.): (61) 

»oo 



xf {t-nu ...,ni) = i-iy-' « [Pn, ■ ... ■ PnM ,l>2. (62) 

Note that there results a correspondence between derivative and variational derivative 
with respect to /„ in TL: 

d S 

Ofn Ofn 

similar to that for P„ (t). Eq.(60) means, roughly speaking, that the l-th order cumulant 
correlator for the gross variables falls off as at — > oo, which expresses the 

expected TL behavior of the thermodynamic fluctuations for the intensive observables. 
It worth emphasizing that adopting Eqs.(40) and Eqs.(60) is by no means sufficient for 
existence of ips {t; f) since it does not guarantee the convergence of the functional Taylor 
series in Eq.(59) 

Finally, consider eh of ip^ {t; f) 

V'(i;f)-lim^,(i;f). (63) 
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Taking formally eh in each term of the series in Eq.(59) with adopting Eq.(40) and using 
the limit 

{t; m) = hm x(^) (t; m) = (P„,)* (64) 
that exists in NESOM-1 unconditionally, gives 

V'(i;f)=V'(^) + E/! J - J ><^'Ht;nu...,ni)l[U,dpi^^, (65) 



provided that the quotients of terms with I > 2 

x«(t;ni,...,nO = lim x« (t; m, = (-!)'"' [P^ 



[P„, • ... • P„J* = hm W [P„, • ... • 

ej.0 



(66) 
(67) 



at least exist. Again, this condition is not sufficient for ip {t; f) to exist. To proceed with 
the proof at goal, the following ultimate statement is postulated. 

Axiom 4 ip {t; f) exists for a set of macro-parameters {including the both BE solutions at 
least) at non-trivial f's in a vicinity of the point f — and is twice functionally differen- 
tiable at that point. 

Remark 1 Thus, the functional ^ (t; f) is not required to be analytic at f = 0. It follows 
with necessity from Axiom 4 and Eqs.(64)-(67) that 



S^l^{t;f) 



Sfn 

S'^ {t; f) 



Sfn^fn 



f=0 



f=0 



{Pnf 
- [Pn ■ Pj 



(68) 
(69) 



but existence of the higher correlators is not necessary for the forthcoming proof. 
4.2 Proof 

Returning now to Eq.(25), put A = {t; f) and B — S {to, to — t) with an arbitrary io- 
Then two inequalities result. The first is 



rj^ J g-5e(t;f) 



S{t,0) -S{to,to-t) 
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Multiplication this inequality by £e^(*° > with subsequent integration over to from 
— oo to i gives 

%m>j:'-^u (70) 

The second inequality is 



Multiplication of this inequality by rje^^^° *\ where 77 > and is arbitrary in any other 
respect, with subsequent integration over from — 00 to i gives 

-00 

t 

+ 6 j Y.Fn{to) [^\Pn{t^.to-t%-^^\Pr:)'°]e<''^-'^dU. (71) 

n 

— oo 

From this point forward, the notation of quasi-invariant is applied to usual functions of 
time, for brevity. 

Dividing by perform TL in Eqs.(70) and (71) with the use of Axioms 1-3, Eqs.(33) 
- (35) and the above stated correspondence between derivatives. This results in 

He it; f) 



Sfn 



(72) 



which does reconfirm non-negativity of [t) = ip^ [t; 0) obtained above, and in 



{t;f) < J ^'HPn)l, fndf^n+ <J, (t) - a (t) 



+ e 



I 



,e{to-t) 



Fn{to) lim (2)(p^(to,to-t))* - (2)(p^)*o diiAdto. 
»oo L ' J I 



When keeping in the latter inequality e finite, t] may freely be tended to zero in its rhs 
without affecting the Ihs which does not depend on r] at all. Let r] to take values of any 
infinitesimal sequence, for which the sequence a {t) converge to its upper limit. This limit 
process, on account of Eq.(36) and Eq.(37), leads to the interim bound 



(i;f) < / P)(P„)7„d//„+ (2)a(i) - lima(i) . 
J nio 



(73) 
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Let e takes values of any infinitesimal sequence, for which the sequence {t) converges 
to its upper limit. Then Eq.(73) results in 



(t;f)= li^V'.(^;f) < / ^^HPnf fndtin + (2)^ (i) - Ih^a^f) (74) 



Let now F„ (io) = ^^^^^(io). Then, in the rhs of Eq(74) (2)(P„)* = (P„)*, while the 
second and third terms cancel each other to give 

V^(i;f)< j {Pnffndlln (75) 

At f = Eq.(75) shows that 

in^^,(0=V^(t;0)<0, 

while this limit should be non negative due to V'e {t) > 0. It is thus proved that on the 
BE solutions of NESOM-2 lim'^^ {t) = 0. So in this case ip does exist and equals zero. It 
is quite remarkable as existence of ip{t;f) has not yet been invoked, but for finishing the 
proof it cannot be avoided. When adopting Axiom 4, Eq.(75) acquires the form 

X{t;f)^ V(*;f) - / {PnYfndfXn< 0. (76) 

and X (^; 0) = V' = 0. This means that, while BE of NESOM-2 are satisfied, the functional 

X {t; f) becomes non-positive attaining atf = absolute maximum {equal zero). Then using 
the necessary condition of functional extremum and Eq.(68) gives 

(t;f) 



Q ^ SX (t; f) 



f=0 



Sfn 



- {p^y = (pj* - (p„)* , (77) 

f=0 



that is BE of NESOM-1. Thus every BE solution of NESOM-2 proves to satisfy also BE 
of NESOM-1, i.e. ^^^F^ (t) — ^^^Fn (t), but not vice versa. That the above extremum is 
indeed maximum is guaranteed by the second variation of % (t; f) at the extremum, i.e. 

S'x* {t; f ) 



5(2)x(t;Af) = J J 



SfJfmdfXndfXn, = - [P (5f) • V (5f)]* < 0, 
f=0 
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where Eq.(69) was used, and V (5f) = J PnSfndjin- The negativity of {t'l ^f) holds 
due to unconditional positivity of the Kubo-Duhamel auto-correlator [8], [13, Preface]. 
This proves the equivalence of NESOM-1 and NESOM-2, if the BE solutions in both 
methods are unique. ■ 

The equivalence in the sense of Eq.(24) is also proven using the techique developed 
above, provided that A e {Ap}, where Ap are some intensive variables (other than gross 
ones), which index p satisfies Axiom 1, and the functional 



■0A (t', q) — — lim lim In Tr 



-5(t,o) -ZgpAp 
e f 



on a set g = {gp} may be constructed to satisfy Axiom 4. In the present framework more 
definite description of variables, for which Eq.(24) holds, seems hardly possible. 

5 Conclusion 

To conclude, the non-pertutbative proof of equivalence between two non-equilibrium en- 
sembles in NESOM, based on MaxEnt principle, is proposed. The proof is thought as 
an improvement upon previous ones [5], [6] and [7, Preface]. Because the rationales of 
NESOM was not clearly delineated, some natural assumptions concerning TL and, inher- 
ent to NESOM, eL are introduced. The present proof, is not also free of some 'extra' 
assumption, namely existence and second-order differentiability of the generating func- 
tional %l) [t] f) However, the latter seems natural in the field-theoretical context and much 
more appropriate than the 'time correlation weackening' [5], [6], entropy-production se- 
ries convergence [6] and 'asymptotic trace normalisation' [7, Preface] conditions, which 
conflict with the basics of NESOM. 

Appendix I 

Consider 'pre-limit' differential form of BE in NESOM-1. Multiplying both sides of 



19 



Eq.(15) by P„ and taking the trace gives 



d_ 

dt 



{Pn)l-^'^(^Pn(t,0)'^ =an,eW, (78) 



where 



(1) 



{Pn-S{t,0))l - « P„- ,5(^,0) 



t 

e 



J [Pn- S {to, to - i))* e'^'°-'^dto. (79) 



-oo 



Thus, for Eq.(22) to hold at a = 1, it is necessary that Zn,e (t) — 0. However, any 
connection of the resulting from this projective-type equations with Eq.(21) can hardly 
be stated. Note that weaker condition 

lim hm Zn,e it) = (80) 

elO |f2|— >oo 

appears in Ref.[5] as a condition for the equivalence of NESOM-1 and NESOM-2. With 
Eq.(80), Eq.(20) for a = 1 holds trivially, but connection of Eq.(19) and that condition 
remains obscure. The key problem here, as seen from Eq.(79), is the time behaviour of ei- 
ther (P„ • Pjn {to, to — t))l or ^Pn- S {to, to — t)^ , which should be at least bounded 
or fall oS at to ^ —oo, respectively. These conditions are not found among the rationales 
of NESOM. 

Appendix II 

This Appendix overviews previous treatments of the equivalence problem. 
Consider first the proof of Ref.[5] . The 'time correlation- weakening' condition of Ref. [5] 
is essentially the same as Eq.(80) and implies that lim lim (Pn- S {to,to — t)) — > 

slO |f!|^oo V /e 

as to —>■ — oo. However, the NESOM perturbational practice evidences that such a behav- 
ior shows up only after BE were used to eliminate the terms oc Fm {to) off S {to, to — t). 
Thus, the condition of Ref. [5] may be necessary, but by no means sufficient, for prooving 
the equivalence. ■ 
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Consider next the proof of Ref.[6] (also [8, Preface]) claimed for classical case. Here 
the approach of Ref.[6] is extended to quantum case. Let us obtain formal expansion of 
^^^Ve i't) with respect the entropy production. Making use of the last equality in Eq.(8), 
expand the operator exponent defining QNESO-1 and its trace to obtain 



(1) 



Pe{t) = (81) 
-oo -oo 



where Vt (A, r) = e^^(*'°) S {t + r, r) e'^^^*'"). Using the identity: S {t + tq, tq) = S {t, 0) - 
/ S (t + Ti,ri) dri, analogous expansion for QNESO-2 is obtained 

TO 



(2) 



Pe{t) = P,(i){l+ (82) 

oo " 1 1 k 

J e'"'^^o y ...y y ...y [Vt{Xun)...-v,{x,,n)ZYidx,dTj 

-oo TO TO 

1 1 

As / ... / [Vt (Ai, Ti) ... • Vt (Afc, rk)]f„ Ilj=i dXj is symmetric in variables n, r^, the in- 



Tl Tk-1 

tegration over them in Eqs.(81) and (82) may be changed to J dri J dT2... J dTk and 

— oo — oo — oo 



/ dTk J dTk-1--- J dri. respectively, to times k\. Using the latter transformation, each term 

TO T/s T2 

of expansion in Eq.(82) can be integrated over tq by parts, after which there remains 
integration over the domain — oo < < Tk-i < ... < ti < 0. Changing this integration 
order to the opposite one makes the integration domains in both Eq.(81) and Eq.(82) the 
same. This results in the unified expansion of QNESO 

^"Ve(t) = P.(t){l + 

OO /*/*/*/*/* ^1 

Ey y - y ^^^^^ y -y [^t(Ai,n)...-v^t(A,,r,)]> h^-^a^' 

^^=^-00 -oo -oo ) 

from which the pre eL quasi- averages expansion follows 

= {A)'+ (83) 

Tl Tfc_l 1 1 



E/ /••• / ^''^9k J... I {[Vt{Xi,n)...-Vt{Xk,n)]l,AyfldX,dT„ 

-oo -oo -oo 
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where '^^^Qk — e^('^i+ - +'^'=\ ^"^^Qk — 6^"^*=. As lim ^"'^Qk — 1, it was concluded in Ref.[6] (also 
[8, Preface]) that Eq.(24) holds unconditionally. In fact, this statement is vahd only if 
the irreducible correlators ([Vi (Ai,ri) ... • Vt (Ajt, rfe)]?^j.yl)* are absolutely integable in the 
domain 0>Ti > ... > Tk > — oo, otherwise counterexamples show that the above integrals 
with a — 1 and a — 2 many have quite different values. This is another form of the 
'time correlation-weakening' condition. Again, referring to the practice of NESOM shows 
that the irreducible correlators, though satisfying a 'spatial correlation-weakening' at all 
Fn{to), don't satisfy 'time correlation- weakening' unless the above-mentioned exclusion 
of 'secular' terms is made with the use of BE. Thus, without the assumption that BE are 
satisfied in advance, the proof of Ref.[6] is deceptive, while that assumption makes the 
proof unclosed. ■ 

At last, consider proof of Ref.[7, Preface]. The use is made of Jensen convexity in- 
equality that gives 

g-5(t,0) > g-5(t,0) 

at each phase-space point. This inequality is consistent with (t) > 0, being a particular 
case of Eq.(70) for f = 0. As a result, in classical case the following pointwise inequality 
between QNESO-2 and QNESO-1 is stated 

^\ (t) - ^\ (t) e-*^W > 0. (84) 

Using the identity 

^'\A)l - ('\A)l ^ J it) - (^)p, it) e-*^W] Adu; + [e'^^^*) - l] ^'\A)1 , 

where dcu is the phase-space measure, and assuming the observable A to be bounded, one 
obtains the estimation 

|(2)(^)* _ < j |(2)p^(i)- Wp,(i)e-*^W||>l|da;+[l-e-*^W] 
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where = max|74|, and (t) > is taken into account. Next, the use of Eq.(84) 

and the trace normahsation of QNESO gives 

W(.l):|<2PL[l-e-*^W] (85) 

Note that this consideration cannot be extended to quantum system, since operator expo- 
nent is not convex operator function [10]. Proceeding with the above estimate, Bitensky 
[7, Preface] assumed two conditions: (i) independence of \\A\\q^ on fl; (ii) the 'asymptotic 
normahsation condition' 

hm hm (t) = 0, or hm (t) = 0. (86) 

elO |f2|— »oo elO 

\n\-*oo 

As seen from Eq.(85), (i) and (ii) lead to the equivalence of the NESOM ensembles in the 
sense of Eq.(24). Although (i) is by no means a restrictive condition, (ii) is in doubt. The 
latter condition is obviously invalid at general F„ (to), because TL and eh of (t) / 
(performed either sequentially or simultaneously) is strictly positive, and so (t) ~^ 
unless BE are satisfied. For TL and sL perfomed sequentially and F„ (to) = ^^^F^ {tg), 
(t) ^ oo, since ip^ (t) remains positive at finite e. Simaltaneous TL and eh of (t) is 
uncertain, since that for (t) / \Q\ is zero. However, for non-stationary process Eq.(86) 
also proves invalid. Indeed, assume that ^^^s exists. Then, given small d > 0, large A > 
and arbitrary r > 0, there exist ^-nd Eq satisfying |Jlo| > ^ (see footnote 1), such 
that 

{to, 0))l - (1)6, (to) > M [a (to) - ('^5 -S]>0 
for every Q D Qq, e < Eq satisfying \Q\ £ > A, and to E [t — r, t\. Making use of Eq.(43) 



and the above inequality gives 

*e (t) > At 



min (7 (to) — 

toe[t-T,t] 



so that ^£ (t) ^ oo as well. For a stationary state, this reasoning is useless, since a — ^^^s 
at Fn — (i)F„. In this case simultaneous TL and eh of remains uncertain, but it 
cannot be arbitrarily fixed at zero value. ■ 
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